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ABSTRACT 

In this work, numerical simulations were used to investigate the gravitational 
stochastic background produced by coalescences occurring up to z ~ 5 of dou¬ 
ble neutron star systems. The cosmic coalescence rate was derived from Monte 
Carlo methods using the probability distributions for forming a massive binary 
and to occur a coalescence in a given redshift. A truly continuous background is 
produced by events located only beyond the critical redshift z* = 0.23. Events 
occurring in the redshift interval 0.027 < z < 0.23 give origin to a “popcorn” 
noise, while those arising closer than z = 0.027 produce a shot noise. The gravita¬ 
tional density parameter Q gw for the continuous background reaches a maximum 
around 670 Hz with an amplitude of 1.1 x 10“ 9 , while the “popcorn” noise has an 
amplitude about one order of magnitude higher and the maximum occurs around 
a frequency of 1.2 kHz. The signal is below the sensitivity of the first generation 
of detectors but could be detectable by the future generation of ground based in¬ 
terferometers. Correlating two coincident advanced-LIGO detectors or two EGO 
interferometers, the expected S/N ratio are respectively 0.5 and 10. 

Subject headings: gravitational waves; neutron star binaries; stochastic back¬ 
ground; laser interferometers 


1. Introduction 


The merger of two neutron stars, two black holes or a black hole and a neutron star 
are among the most important sources of gravitational waves (GW), due to the huge energy 
released in the process. In particular, the coalescence of double neutron stars (DNS) may 
radiate about 10 53 erg in the last seconds of their inspiral trajectory, at frequencies up to 
1.4-1.6 kHz, range covered by most of the ground-based laser interferometers like VIRGO 
(Bradaschia et al. 1990), LIGO (Abramovici et al. 1992), GEO (Hough 1992) or TAMA 
(Kuroda et al. 1997). Besides the amount of energy involved in these events, the rate at which 
they occur in the local universe is another parameter characterizing if these mergings are or 
not potential interesting sources of GW. In spite of the large amount of work performed in 
the past years, uncertainties persist in estimates of the DNS coalescence rate. In a previous 
investigation, we have revisited this question (de Freitas Pacheco et al. 2005; Regimbau 
et al. 2005), taking into account the galactic star formation history derived directly from 
observations and including the contribution of elliptical galaxies when estimating the mean 
merging rate in the local universe. Based on these results, we have predicted a detection 
rate of one event every 125 and 148 years by initial LIGO and VIRGO respectively and up 
to 6 detections per year in their advanced configurations. 

Besides the emission produced by the coalescence of the nearest DNS, the superposition 
of a large number of unresolved sources at high redshifts will produce a stochastic background 
of GW. In the past years, different astrophysical processes susceptible to generate a stochastic 
background have been investigated. On the one hand, distorted black holes (Ferrari et 
al. 1999; de Araujo et al. 2000), bar mode emission from young neutron stars (Regimbau 
2001) are examples of sources able to generate a shot noise (time interval between events 
large in comparison with duration of a single event), while supernovas or hypernovas (Blair 
& Ju 1997; Coward et al. 2001, 2002; Buonanno et al. 2004) are expected to produce an 
intermediate “popcorn” noise. On the other hand, the contribution of tri-axial rotating 
neutron stars (Regimbau & de Freitas Pacheco 2001), including magnetars (Regimbau & de 
Freitas Pacheco 2005), constitutes a truly continuous background. 

Populations of compact binaries as, for instance, the cataclismic variables are respon¬ 
sible for the existence of a galactic background of GW in the mHz domain, which could 
represent an important source of confusion noise for space detectors as LISA (Evans et al. 
1987; Hils et al. 1990; Bender & Hils 1997; Postnov & Prokhorov 1998; Nelcnan et al. 2001; 
Timpano et al. 2005). These investigations have been extended recently to the extra-galactic 
contribution. Schneider et al. (2001), Kim et al. (2003) and Cooray (2004) considered cos¬ 
mological populations of double and mixed systems involving black holes, neutron stars and 
white dwarfs, while close binaries originated from low and intermediate mass stars were 
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discussed by Farmer & Phinney (2002). 

In this work, using the DNS merging rate estimated in our precedent study, we have 
estimated the gravitational wave background spectrum produced by these coalescences. Nu¬ 
merical simulations based on Monte Carlo methods were performed in order to determine 
the critical redshift z * beyond which the duty cycle condition required to have a continuous 
background (D > 1) is satisfied. Unlike previous studies which focus their attention on the 
early low frequency inspiral phase covered by LISA (Schneider et al. 2001; Farmer & Phinney 
2002; Cooray 2004), here we are mainly interested in the few thousand seconds before the 
last stable orbit is reached, when more than 96% of the gravitational wave energy is released. 
The signal frequency is in the range 10-1500 Hz, covered by ground based interferometers. 
The paper is organized as follows. In §2, the simulations are described; in §3 the contribution 
of DNS coalescences to the stochastic background is calculated; in §4 the detection possibil¬ 
ity with laser beam interferometers is discussed and, finally, in §5 the main conclusions are 
summarized. 


2. The Simulations 


In order to simulate by Monte Carlo methods the occurrence of merging events, we have 
adopted the following procedure. The first step was to estimate the probability for a given 
pair of massive stars, supposed to be the progenitors of DNS, be formed at a given redshift. 
This probability distribution is essentially given by the cosmic star formation rate (Coward 
et al. 2002), normalized in the redshift interval 0 < z < 5, e.g., 


Pi w = 


dRf(z)/dz 


N n 


( 1 ) 


The normalization factor in the denominator is essentially the rate at which massive binaries 
are formed in the considered redshift interval, e.g., 


N p — ( dRf(z) / dz)dz , 


( 2 ) 


which depends on the adopted cosmic star formation rate, as we shall see later. 


The formation rate of massive binaries per redshift interval is 


Rzf( z f) 


dR,(z,) , R}(z f ) dV(z f ) 

- — - 

dzf 1 + Zf dzf 


(3) 


In the equation above, R*f(z) is the cosmic star formation rate (SFR) expressed in M 0 
Mpc _3 yr _1 and X p is the mass fraction converted into DNS progenitors. Hereafter, rates per 



4 


comoving volume will always be indicated by the superscript while rates with indexes 
“z f v or “z c ” refer to differential rates per redshift interval, including all cosmological factors. 
The (1+z) term in the denominator of eq. 3 corrects the star formation rate by time dilatation 
due to the cosmic expansion. In the present work we assume that the parameter X p does not 
change significantly with the redshift and thus it will be considered as a constant. In fact, 
this term is the product of three other parameters, namely, 

Ap = ^Nsfb^NS (4) 

where (5ns is the fraction of binaries which remains bounded after the second supernova 
event, fb is the fraction of massive binaries formed among all stars and Xns is the mass 
fraction of neutron star progenitors. 

According to the results by de Freitas Pacheco et al. (2005),Regimbau et al. (2005), 
/ 3ns = 0.024 and fb = 0.136, values which will be adopted in our calculations. Assuming 
that progenitors with initial masses above 40 M 0 will produce black holes and considering 
an initial mass function (IMF) of the form £(m) = dm -7 , with 7 = 2.35 (Salpeter’s law), 
normalized within the mass interval 0.1 - 80 M 0 such as f m£(m)dm = 1, it results finally 
A ns = fg° €(m)dm = 5.72 x 10 ~ 3 Mq 1 and X p = 1.85 x 10 ~ 5 Mg 1 . The evaluation of the 
parameters (5ns and fb depends on different assumptions, which explain why estimates of 
the coalescence rate of DNS found in the literature may vary by one or even two orders of 
magnitude. The evolutionary scenario of massive binaries considered in our calculations (see 
de Freitas Pacheco et al. (2005) for details) is similar to that developed by Bclczynski et 
al. ( 2002 ), in which none of the stars ever had the chance of being recycled by accretion. 
Besides the evolutionary path, the resulting fraction (5ns of bound NS-NS binaries depends 
on the adopted velocity distribution of the natal kick. The imparted kick may unbind binaries 
which otherwise might have remained bound or, less probably, conserve bound systems which 
without the kick would have been disrupted. The adopted value for (5ns corresponds to a 
1-D velocity dispersion of about 80 km/s. This value is smaller than those usually assumed 
for single pulsars, but consistent with recent analyses of the spin period-eccentricity relation 
for NS-NS binaries (Dewi et al. 2005). Had we adopted a higher velocity dispersion in our 
simulations (230 km/s instead of 80 km/s), the resulting fraction of bound systems is reduced 
by one order of magnitude, e.g., (5ns = 0.0029. If, on the one hand the fraction of bound NS- 
NS systems after the second supernova event depends on the previous evolutionary history 
of the progenitors and on the kick velocity distribution, on the other hand estimates of the 
fraction fb of massive binaries formed among all stars depends on the ratio between single 
and double NS systems in the Galaxy and on the value of (5ns itself (de Freitas Pacheco et al. 
2005). We have estimated relative uncertainties of about <jf b / fb ~ 0.5 and vp NS /(5 ns ~ 0.75, 
leading to a relative uncertainty in the parameter X p of about cr\ p /X p ~ 0.9. However, we 
emphasize that these are only formal uncertainties resulting from our simulations, which 
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depend on the adopted evolutionary scenario for the progenitors. A comparison with other 
estimates can be found in de Freitas Pacheco et al. (2005). 

The element of comoving volume is given by 


dV (z) = 47rr(z) 2 


c dz 
lf 0 E(Sli,z) 


(5) 


with 

E{n u z) = [fi m (l + ^) 3 + ^] 1/2 (6) 

where and £l v are respectively the present values of the density parameters due to matter 
(baryonic and non-baryonic) and vacuum, corresponding to a non-zero cosmological constant. 
A “flat” cosmological model + = 1) was assumed. In our calculations, we have taken 

Q m = 0.30 and Q„ = 0.70, corresponding to the so-called “concordance” model derived from 
observations of distant type la supernovae (Schmidt et al. 1998) and the power spectra of 
the cosmic microwave background fluctuations (Spergel 2003). The Hubble parameter H 0 
was taken to be 65 kms -1 Mpc -1 . 

Porciani & Madau (2001) provide three models for the cosmic SFR history up to redshifts 
z ~ 5. Differences among these models are mainly due to various corrections applied, 
in particular those due to extinction by the cosmic dust. In our computations, we have 
considered the second model, labelled SFR2 (Madau et al. 1998) but numerical results using 
SFR1 (Steidel 1999) will be also given for comparison. Both rates increase rapidly between 
z ~ 0 — 1 and peak at z ~ 1 — 2, but SFR1 decreases gently after z ~ 2 while SFR2 remains 
more or less constant (Figure 1). 

The next step consists to estimate the redshift z b at which the progenitors have already 
evolved and the system is now constituted by two neutron stars. This moment fixes also 
the beginning of the inspiral phase. If p, (« 10 8 yr) is the mean lifetime of the progenitors 
(average weighted by the IMF in the interval 9-40 M 0 ) then 


Zb = z f - H 0 r b ( 1 + z f )E(z f ) (7) 

Once the beginning of the inspiral phase is established, the redshift at which the coa¬ 
lescence occurs is estimated by the following procedure. The duration of the inspiral phase 
depends on the orbital parameters just after the second supernova and on the neutron star 
masses. The probability for a given DNS system to coalesce in a timescale r was initially 
derived by de Freitas Pacheco (1997), confirmed by subsequent simulations (Vincent 2002; 
de Freitas Pacheco et al. 2005) and is given by 



6 


P t {t) = B/t 


( 8 ) 


Simulations indicate a minimum coalescence timescale To = 2 x 10 5 yr but a considerable 
number of systems have a coalescence timescale higher than the Hubble time. The normalized 
probability in the range 2 x 10 5 yr up to 20 Gyr implies B = 0.087. Therefore, the redshift 
z c at which the coalescence occurs after a timescale r is derived from the equation 


HoT = L (1 + 4B(i) 


( 9 ) 


which was solved in our code by an iterative method. The resulting distribution of the 
number of coalescences as a function of z c is shown in Figure 2 for both star formation rates 
SFR1 and SFR2, while the corresponding coalescence rate per redshift interval, R zc (z), is 
shown in Figure 1. In the same figure, for comparison, we have plotted the formation rate 
R z f(z) (eq. 3). Notice that the maximum of R zc (z ) is shifted towards lower redshifts with 
respect to the maximum of R z f{z ), reflecting the time delay between the formation of the 
progenitors and the coalescence event. The coalescence rate R zc (z), does not fall to zero at 
z — 0 because a non negligible fraction of coalescences (~ 3% for SFR2 and ~ 5% for SFR1) 
occurs later than z — 0. 


3. The Gravitational Wave Background 


The nature of the background is determined by the duty cycle defined as the ratio of the 
typical duration of a single burst f to the average time interval between successive events, 


e-g., 


D(z*) = r(l + z')R zc (z')dz' 


( 10 ) 


Jo 

The critical redshift z* at which the background becomes continuous is fixed by the condition 
D(z) > 1. Since we are interested in the last instants of the inspiral, when the signal is 
within the frequency band of ground based interferometers, we took f = 1000 s, duration 
which includes about 96% of the total energy released (see Table 1). From our numerical 
experiments and imposing D = 1, one obtains z* = 0.23 when the SFR2 is used and z* = 0.27 
for the SFR1. About 96% (94% in the case of SFR1) of coalescences occur above such a 
redshift, contributing to the production of a continuous background. Sources in the redshift 
interval 0.027 < z < 0.23 (SFR2) or 0.032 < z < 0.27 (SFR1) correspond to a duty cycle D 
= 0.1 and they are responsible for a cosmic “popcorn” noise. 


The gravitational fluence (given here in erg cm 2 Hz x ) in the observer frame produced 
by a given DNS coalescence is: 
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Uo = 


1 (IE, 


gw 


4:ird 2 L 


du 


(1 + z c ) 


( 11 ) 


where cU = (1 + z c )r is the distance luminosity, r is the proper distance, which depends on 
the adopted cosmology, dE gw /dv is the gravitational spectral energy and v = (1 + z c )v Q the 
frequency in the source frame. In the quadrupolar approximation and for a binary system 
with masses rni and m 2 in a circular orbit: 


dE gw /du = Kv-' l% 


( 12 ) 


where the fact that the gravitational wave frequency is twice the orbital frequency was taken 
into account. Then 

(Gbr) 2 / 3 mim 2 

3 (mi + m 2 ) 1 / 3 

Assuming rn\ = m 2 = 1.4 one obtains K = 5.2 x 10 50 erg Hz -2 / 3 . 

The spectral properties of the stochastic background are characterized by the dimen¬ 
sionless parameter (Ferrari et al. 1999): 



Qgwiyo) 3 


C Pc 




(14) 


where v D is the wave frequency in the observer frame, p c is the critical mass density needed 
to close the Universe, related to the Hubble parameter H 0 by, 


Pc = 


3 H 2 0 
8tt G 


(15) 


F Uo is the gravitational wave flux (given here in erg cm 2 Hz x s x ) at the observer frequency 
u a , integrated over all sources at redshifts z c > z*, namely 


Fu a = 




(16) 


Instead of solving analytically the equation above by introducing, for instance, an ade¬ 
quate fit of the cosmic coalescence rate, we have calculated the integrated gravitational flux 
by summing individual fluences (eq. 11), scaled by the ratio between the expected number 
of events per unit of time and the number of simulated coalescences or, in other words, the 
ratio between the total formation rate of progenitors (eq. 2) and the number of simulated 
massive binaries, e.g., 



N,, 


j vr iy svm 

F v „ = p 


N. 




i= 1 


( 17 ) 


The number of runs (or N S i m ) in our simulations was equal to 10 7 , representing an 
uncertainty of < 0.1% in the density parameter Q gw . Using the SFR2, the derived formation 
rate of progenitors is N p = 0.031 s^ 1 , whereas for the SFR1 one obtains N p = 0.024 s _1 . For 
each run, the probability distribution Pf(z) defines, via Monte Carlo, the redshift at which 
the massive binary is formed. The beginning of the inspiral phase at Zf, is fixed by eq. 7. 
Then, in next step, the probability distribution of the coalescence timescale and eq. 9 define 
the redshift z c at which the merging occurs. The fluence produced by this event is calculated 
by eq. 11, stored in different frequency bins in the observer frame and added according to 
the equation above. 

Figure 4 shows the density parameter Q gw as a function of the observed frequency 
derived from our simulations. The density parameter il gw increases as v^J 3 at low frequencies, 
reaches a maximum amplitude of about 1.1 x 1CT 9 around 670 Hz in the case of SFR2 and 
a maximum of 8.4 x 10 _1 °) around 630 Hz, in the case of SFR1. A high frequency cut-off at 
~ 1220 Hz (~ 1170 Hz for SFR1 ) is observed, corresponding approximately to the frequency 
of the last stable orbit at the critical redshift 2* = 0.23 (2* = 0.27 for SFR1). Calculations 
performed by Schneider et al. (2001), in spite of the similar local merging rates, indicate 
that the maximum occurs at lower frequencies (~ 100 Hz) with an amplitude (scaled to the 
Hubble parameter adopted in this work) lower by a factor of seven. However, as those authors 
have stressed, their calculations are expected to be accurate in the frequency range 10 /xHz 
up to 1 Hz, since they have set the value of the maximum frequency v max to about that 
expected at a separation of three times the “last stable orbit” (LSO), e.g., u max ~ 0.19 vlso- 
Thus, a direct comparison with our results is probably meaningless. 

A more conservative estimate can be obtained if one adopts a higher duty cycle value, 
namely, D > 10, corresponding to sources located beyond 2 ~ 1.05 (2 ~ 1.075 for SFR1). 
Our results are plotted in Fig. 4, which includes, for comparison, the “popcorn” noise con¬ 
tribution arising from sources between 0.027 < 2 < 0.23 (0.032 < 2 < 0.27 for SFR1), 
corresponding to the intermediate zone between a shot noise (D < 0.1) and a continuous 
background (D > 1). When increasing the critical redshift (or removing the nearest sources), 
the amplitude of Q gw decreases and the spectrum is shifted toward lower frequencies. The 
amplitude of the “popcorn” background is about one order of magnitude higher than the 
continuous background, with a maximum of about Q gw = 1.3 x 1(T 8 (8.8 x 10" 9 for SFR1) 
around a frequency of 1.2 kHz. 

Since some authors use, instead of Q gw , the gravitational strain ydefined by Allen 
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& Romano (1999) as 

S„K) = ( 18 ) 

we show this quantity in Figure 5. 


4. Detection 

Because the background obeys a Gaussian statistic and can be confounded with the 
instrumental noise background of a single detector, the optimal detection strategy is to 
cross-correlate the output of two (or more) detectors, assumed to have independent spectral 
noises. The cross correlation product is given by (Allen & Romano 1999): 


/ OO 

8~1 *(f)Q(m(f)df (19) 

-OO 

where 

Q(f) oc r (f) n 9w(f) (20) 

is a filter that maximizes the signal to noise ratio ( S/R ). In the above equations, Pi(f) and 
P 2 (_/) are the power spectral noise densities of the two detectors and T is the non-normalized 
overlap reduction function, characterizing the loss of sensitivity due to the separation and 
the relative orientation of the detectors. The optimized S/N ratio for an integration time T 
is given by (Allen 1997) : 


= sh[ t r r r VHM) 

n 1 8 tt 4 j 0 1 pp^mur 


( 21 ) 


In the literature, the sensitivity of detector pairs is usually given in terms of the minimum 
detectable amplitude for a flat spectrum (tt gw equal to constant) (Allen & Romano 1999), 

e-g-, 


At ,- 2 r°° r 2 (f) 

Qmin = 3 Hf/T (eriC ~ 1{2a) " erfc “ 1(27))[ / df pP 1 (f)P 2 (f) ] ~ 1/2 (22) 

The expected minimum detectable amplitude for the main pair of detectors in the world, 
after one year integration, are given in Table 2, for a detection rate a = 90% and a false alarm 
rate 7 = 10%. The power spectral densities expressions used for the present calculation 
can be found in Darnour et al. (2001). Q ?mn is of the order of 10 -6 — 10 ~ 5 for the first 
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generation of interferometers combined as LIGO/LIGO and LIGO/VIRGO. Their advanced 
counterparts will permit an increase of two or even three orders of magnitude in sensitivity 
(Q mm ~ 1CT 9 — 10” 8 ). The pair formed by the co-located and co-aligned LIGO Hanford 
detectors, for which the overlap reduction function is equal to one, is potentially one order 
of magnitude more sensitive than the Hanford/Livingston pair, provided that instrumental 
and environmental noises could be removed. 

However, because the spectrum of DNS coalescences is not flat and the maximum oc¬ 
curs out of the optimal frequency band of ground based interferometers, which is typically 
around 50 — 300 Hz, as shown in Figure 6, the S/N ratio is slightly reduced. Considering 
the co-located and co-aligned LIGO interferometer pair, we find a signal-to-noise ratio of 
S/N ~ 0.002 ( S/N ~ 0.5) for the initial (advanced) configuration. Unless the coalescence 
rate be substantially higher than the present expectations, our results indicate that their 
contribution to the gravitational background is out of reach of the first and the second gen¬ 
eration of interferometers. On the other hand, the sensitivity of the future third generation 
of detectors, presently in discussion, could be high enough to gain one order of magnitude 
in the expected S/N ratio. Examples are the Large Scale Cryogenic Gravitational Wave 
Telescope (LCGT), sponsored by the LIniversity of Tokyo and the European antenna EGO 
(Sathyaprakash, private communication). EGO will incorporate signal recycling, diffractive 
optics on silicon mirrors, cryo-techniques and kW-class lasers, among other technological 
improvements. A possible sensitivity for this detector is shown in Figure 5, compared to the 
expected sensitivity of advanced LIGO. Around 650 Hz, the planned strain noise S n (u) 
is about 8 x 10 -24 Hz -1 / 2 for the advanced LIGO configuration while at this frequency, the 
planned strain noise for EGO is 2 x 10 -24 Hz -1 / 2 , which represents a gain by a factor of ~ 4. 
Considering two interferometers located at the same place, we find a signal-to-noise ratio 
S/N ~ 10. 

On the other hand, the popcorn noise contribution could be detected by new data 
analysis techniques currently under investigation, such as the search for anisotropies (Allen 
& Ottewill 1997) that can be used to create a map of the GW background (Cornish 2001), the 
maximum likelihood statistic (Drasco & Flanagan 2003), or methods based on the Probability 
“Event Horizon” concept (Coward & Burman 2005), which describes the evolution, as 
a function of the observation time, of the cumulated signal throughout the LIniverse. The 
PEH of the GW signal evolves fastly from contributions of high redshift populations, forming 
a real continuous stochastic background, to low redshift and less probable sources that can 
be resolved individually, while the PEH of the instrumental noise is expected to evolve much 
slower. Consequently, the GW signature could be distinguished from the instrumental noise 
background. 



5. Conclusions 


In this work, we have performed numerical simulations using Monte Carlo techniques to 
estimate the occurrence of double neutron star coalescences and the gravitational stochastic 
background produced these events. Since the coalescence timescale obeys a well defined 
probability distribution (P(r) oc 1/r), derived from simulations of the evolution of massive 
binaries (de Freitas Pacheco et al. 2005), the cosmic coalescence rate does not follow the 
cosmic star formation rate and presents necessarily a time-lag. In the case where the sources 
are supernovae or black holes, the gravitational burst is produced in a quite short timescale 
after the the formation of the progenitors. Therefore, the time-lag is negligible and the 
comoving volume where the progenitors are formed is practically the same as that where 
the gravitational wave emission occurs, introducing a considerable simplification in the cal¬ 
culations. This is not the case when NS-NS coalescences are considered, since timescales 
comparable or even higher than the Hubble timescale have non negligible probabilities. The 
maximum probability to form a massive binary occurs at z ~ 1.7, depending slightly on the 
adopted cosmic star formation rate, whereas the maximum probability to occur a coalescence 
is around z 

We have found that a truly continuous background is formed only when sources located 
beyond z > 0.23 (z > 0.27 for the SFR1 case), including 96% (94% for SFR1) of all events 
and the critical redshift corresponds to the condition D > 1. Sources in the redshift interval 
0.027 < z < 0.23 (0.032 < z < 0.27 for SFR1) produce a “popcorn” noise. Our computations 
indicate that the density parameter Q gw has a maximum around 670 Hz (630 Hz for SFR1), 
attaining an amplitude of about of 1.1 x 10 9 (8.3 x 10 10 for SFR1). The low frequency 
cutoff around 1.2 kHz corresponds essentially to the gravitational redshifted wave frequency 
associated to last stable orbit of sources located near the maximum of the coalescence rate. 

The computed signal is below the sensitivity of the first and the second generation of 
detectors. However, using the planned sensitivity of third generation interferometers, we 
found that after one year of integration, the cross-correlation of two EGO like coincident 
antennas, gives an the optimized signal-to-noise of S/R ~ 10. 

The “popcorn” contribution is one order of magnitude higher with a maximum of 
f } gw ~ 1.3 x 10~ 8 (8.8 x 10“ 9 for SFR1) at ~ 1.2 kHz. This signal, which is character¬ 
ized by the spatial and temporal evolution of the events as well as by its signature, can be 
distinguished from the instrumental noise background and adequate data analysis strategies 
for its detection are currently under investigation (Allen & Ottewill 1997; Cornish 2001; 
Drasco & Flanagan 2003; Coward & Burman 2005). 
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Table 1. For three values of the emission frequency (column 1), the corresponding time left 
to the last stable orbit (second column) and percentage of energy released (third column) 


Vmin (Hz) 

T 

A E/E t (%) 

100 

2 s 

84 

10 

1000 s 

96 

1 

5.26 d 

99 
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Table 2. Expected il m in of different interferometer pairs for a flat background spectrum 
and an integration time T = 1 year, a detection rate a = 90% and a false alarm rate 
7 = 10%. LHO and LLO stand for LIGO Hanford Observatory and LIGO Livingston 

Observatory 



LHO-LHO 

LHO-LLO 

LLO-VIRGO 

VIRGO-GEO 

initial 

advanced 

4 x 10” 7 

6 x 10~ 9 

4 x 10“ 6 

1 x 10 -9 

8 x 10“ 6 

8 x 10~ 6 



Fig. 1.— Cosmic star formation rate Rf(z ) in M 0 Mpc~ 3 yr _1 in a flat cosmology with Q m 
= 0.30 and 0,, =0.70: SFR2 is represented by a black line and SFR1 by a grey line. The 
formation and the coalescence rate per redshift interval are also plotted for comparison. 




number of coalescence 
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Fig. 2.— Distribution of coalescences as a function of the redshift derived from our simula¬ 
tions, corresponding to 10' numerical experiments. 























































































































19 



Fig. 3.— Spectrum of the expected gravitational energy density parameter Q gw for the 
continuous regime. Results are shown for the two star formation rates adopted in this work. 
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Fig. 4.— Spectrum of the expected gravitational energy density parameter Q gw correspond¬ 
ing to NS-NS coalescences occurring beyond z* = 0.23 (bold continuous curve). A more 
conservative background (D > 10) corresponding to sources beyond z = 1.05 is plotted 
for comparison (grey continuous curve). The expected “popcorn” noise contribution arising 
from sources in the redshift interval 0.027 < z < 0.23 (dashed curve) is also shown. This 
corresponds to the intermediate zone between a shot noise (D < 0.1) and a continuous 
background (D > 1). 
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Fig. 5.— Gravitational strain in Hz -1 / 2 corresponding to NS-NS coalescences occurring 
beyond z* = 0.23, along with the planned sensitivity curves of LIGO Ad (continuous curve) 
and EGO (dashed curve). 
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Fig. 6.— Integrand of the signal-to-noise ratio (eq. 21) for advanced LIGO pairs: LHO-LHO 
(grey curve) and LHO-LLO (black curve) 



















